function [clusterIdx, H_normalized] = Coreg_config_auto(Ks, nCluster)
[nSmp, ~, nKernel] = size(Ks);

options = [];
options.lambda=0.01;
H = coreg_centroid_on_multi_kernel(Ks, nCluster, options);
H_normalized = H ./ repmat(sqrt(sum(H.^2, 2)), 1,nCluster);
clusterIdx = litekmeans(H_normalized, nCluster, 'MaxIter',100, 'Replicates', 50);